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Abstract.  Lunar  occupations  were  analyzed  to  detect  errors  in  the 
proper  motion  system  of  the  Hipparcos  frame.  The  analysis  suggests 
much  larger  errors  in  the  rotation  of  the  Hipparcos  reference  frame  than 
those  explained  by  the  Hipparcos  team,  which  were  determined  from  the 
direct  comparison  between  the  Hipparcos  and  FK5  catalogs. 


1.  Introduction 

The  Fundamental  Catalog  FK5  (Fricke  et  al.,  1988)  had  been  used  for  positions 
of  objects  on  the  sky,  but  it  is  well-known  that  the  IAU  1976  precession  constant 
used  in  the  FK5  has  an  error  A p  of  about  —3.0  mas/year,  which  had  been 
independently  obtained  from  VLBI  (McCarthy  &  Luzum,  1991),  lunar  laser 
ranging  (Williams  et  al.,  1991),  and  proper  motion  analyses  (Miyamoto  &  Soma, 
1993).  Therefore  one  can  expect  that  the  FK5  system  should  have  the  following 
rotation  with  respect  to  the  ICRS  (International  Celestial  Reference  System): 

uj x  ~  0 

ujy  =  —  Apsine  ~  +1.2  mas/year  (1) 

ojz  =  Apcose  —  AE  ~  —1.6  mas/year 

where  we  use  the  value  of  the  equinox  motion  A E  of  about  —1.2  mas/year  by 
Miyamoto  &  Soma  (1993). 

The  Hipparcos  Catalogue  (ESA,  1997)  was  released  in  1997.  The  Hippar¬ 
cos  reference  frame  is  said  to  be  linked  to  the  ICRS  with  the  accuracy  of  0.6 
mas  (milliarcseconds)  in  orientation  and  0.25  mas/year  in  spin,  but  the  direct 
comparison  of  the  proper  motions  between  Hipparcos  and  FK5  made  by  Feissel 
and  Mignard  (1998)  gives  the  following  rotation  of  the  FK5  system  with  respect 
to  the  Hipparcos  system: 


^ X 

=  -0.30 

mas /year 

LOy 

=  +0.60 

mas /year 

(2) 

A 

=  +0.70 

mas /year 

which  is  inconsistent  with  equation  (1). 

Therefore  it  is  important  to  evaluate  the  accuracy  of  the  Hipparcos  proper 
motion  system  by  other  methods.  Now  that  the  lunar  positions  in  the  latest 
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JPL  lunar  and  planetary  ephemerides  DE405  (Standish,  1998)  have  mas-level 
accuracy  with  respect  to  the  ICRS  because  of  the  very  high  accuracy  of  the 
lunar  laser  ranging  measurements,  the  problem  can  be  resolved  by  measuring  the 
lunar  positions  with  respect  to  the  Hipparcos  stars.  The  author  analyzed  lunar 
occultation  observations  made  from  1955  using  the  DE405  ephemerides  and 
the  Hipparcos  and  ACT  catalogs  to  detect  the  possible  errors  in  the  Hipparcos 
proper  motion  system.  This  paper  discusses  the  results  of  the  analysis. 


2.  Observations  of  Lunar  Occultations 

Lunar  occultation  observations  were  collected  until  1980  by  RGO  (Royal  Green¬ 
wich  Observatory)  and  after  that  by  ILOC  (International  Lunar  Occultation 
Center)  in  Tokyo.  The  ILOC  released  the  observations  made  by  the  end  of 
1997.  Since  it  is  important  to  convert  the  observed  Universal  Time  to  a  uniform 
time  scale  (TAI  in  this  case)  in  the  analysis,  only  the  observations  made  after 
1955.5  are  analyzed  in  this  paper.  In  addition,  visual  observations  suffer  from 
personal  equations  (observers’  reaction  times)  and  can  give  secular  effects  in  the 
analysis.  Therefore  we  only  deal  with  observations  made  photoelectrically  or  by 
video  equipment. 

The  total  number  of  such  observations  is  15,430.  The  number  of  such 
observations  in  each  year  is  given  in  Table  1. 


Table  1.  Number  of  Photoelectric  and  Video  Observations. 


Year 

No. 

Year 

No. 

Year 

No. 

Year 

No. 

Year 

No. 

1961 

60 

1971 

369 

1981 

528 

1991 

549 

1962 

15 

1972 

517 

1982 

571 

1992 

674 

1963 

11 

1973 

516 

1983 

779 

1993 

619 

1964 

33 

1974 

395 

1984 

611 

1994 

543 

1955 

10 

1965 

33 

1975 

515 

1985 

557 

1995 

904 

1956 

57 

1966 

26 

1976 

502 

1986 

471 

1996 

616 

1957 

34 

1967 

63 

1977 

695 

1987 

478 

1997 

633 

1958 

21 

1968 

97 

1978 

559 

1988 

395 

1959 

21 

1969 

200 

1979 

393 

1989 

467 

1960 

59 

1970 

260 

1980 

143 

1990 

431 

3.  Analysis 

The  analysis  were  made  using  the  following  parameters  based  on  the  analysis 
by  Soma  (1985): 

Corrections  to  the  Moon’s  positions:  AA  =  +0"50  and  A/3  =  — 0"24 
Radius  of  the  Moon:  i?o  =  1738.107  km 

Correction  for  the  ellipticity  of  Watts’  datum:  A p  =  +0'/128  sin  2(WA) 

Relation  of  Watts  Angle  (WA)  with  Axis  Angle  (AA):  WA  =  AA  +  0°.24 

The  values  AA  and  A/3  above  are  the  corrections  to  the  Moon’s  ecliptic  longi¬ 
tude  and  latitude  at  mean  distance  (384, 400  km)  in  order  to  compensate  the 
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difference  between  the  center  of  the  reference  datum  of  Watts’  lunar  limb  profile 
charts  and  the  center  of  mass  of  the  Moon. 

The  lunar  limb  profile  irregularity  was  removed  by  the  corrections  given 
by  Watts  (1963).  In  1994,  lunar  topographic  features  were  observed  by  the 
Clementine  lunar  mapping  mission  (Smith  et  al. ,  1997),  but  the  Clementine 
laser  altimeter  data  are  too  sparse  (they  only  covered  narrow  meridional  strips 
separated  by  almost  3°  of  selenographic  longitude)  to  replace  Watts’  data,  and 
the  altimeter  data  for  the  polar  regions  were  not  obtained.  Observations  not 
having  Watts’  limb  profile  corrections  were  rejected  because  these  can  have  sys¬ 
tematic  effects  on  the  solution. 

The  lunar  positions  were  calculated  from  JPL’s  latest  ephemerides  DE405 
(Standish,  1998).  The  positions  of  stars  were  taken  from  the  Hipparcos  Cata¬ 
logue  (ESA,  1997).  For  the  stars  not  given  in  the  Hipparcos  Catalogue  the  ACT 
catalog  (Urban,  1998)  was  used.  The  lunar  and  stellar  positions  thus  calculated 
are  all  claimed  to  be  on  ICRS. 

The  conversion  from  ICRS  to  MCEP  (Earth’s  mean  equator  and  mean 
equinox  of  J2000)  was  made  using  the  parameters  determined  by  Chapront  et 
al.,  (1999).  Precession  from  J2000  and  sidereal  time  were  calculated  using  the 
formulas  by  Williams  (1994)  that  were  based  on  the  general  precession  correction 
of  —3.266  mas/year  to  the  IAU  1976  value.  The  nutation  was  from  the  ZMOA 
1990  solution  by  Herring  (1991).  The  values  of  UT1  —  UTC  and  the  data  for 
the  Earth’s  polar  motion  were  from  BIH  and  IERS.  The  geodetic  coordinates 
based  on  various  local  geodetic  datums  were  transformed  to  WGS  84  using  the 
transformation  parameters  given  by  NIMA  (1997).  The  geoid  heights  based  on 
EGM96  were  taken  into  account. 

The  height  A  a  of  the  star  from  the  lunar  limb,  calculated  from  the  observed 
times  were  used  to  solve  the  correction  AA  to  the  lunar  ecliptic  longitude,  the 
correction  Aa  to  the  right  ascension  system  of  the  stars.  The  rotation  Ac  of  the 
stars’  coordinates  around  the  equinox  including  their  derivatives  with  respect  to 
time  were  obtained  from  the  following  equation: 

Act  =  A  A  cos  (3  ±  Aci  cos  A  sin  e  cos  N  ±  Ae(cos  A  sin  (3  sin  N  —  sin  A  cos  N),  (3) 

where  A  and  (3  are  the  ecliptic  longitude  and  latitude  of  the  Moon,  e  is  the 
obliquity  of  the  ecliptic,  and  N  is  the  position  angle  of  the  event  measured  from 
the  projection  of  the  ecliptic  north  pole. 

The  number  of  equations  is  14,  963  and  the  mean  epoch  is  1983.9.  The  RMS 
of  the  residuals  is  0^201.  The  solution  in  mas  is  as  follows: 

AA  =  -94.7  ±0.56(Y-  1980), 

±  2.2  ±0.22 

Aa  =  ±48.4  -  4.06(Y-  1980),  (4) 

±12.7  ±  1.27 

Ac  =  -17.0  ±0.67(Y-  1980). 

±  4.2  ±0.43 

where  Y  is  the  year.  The  values  at  the  epoch  1991.25  of  the  Hipparcos  Catalogue 
are  found  to  be 

AA  =  —88  ±  3  mas, 
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A  a 
Ae 


+  3  ±  19  mas, 
—  9+6  mas. 


(5) 


The  constant  term  in  AA  can  be  regarded  as  the  correction  to  the  value  of  the 
difference  between  the  center  of  the  reference  datum  of  Watts’  charts  and  the 
center  of  mass  of  the  Moon.  Since  the  value  +0"50  was  adopted  for  it,  the 
solution  suggests  that  the  longitude  component  of  the  difference  of  the  centers 
is  +0^41  in  the  sense  that  the  center  in  Watts’  charts  leads  the  center  of  mass 
in  the  lunar  orbit.  The  constant  terms  in  Aa  and  Ae  are  found  to  be  not  very 
significant. 


The  linear  terms  suggest  the  following 
frame: 

rotation  of  the 

Hipparcos 

reference 

^ X 

=  (Ae)t 

+0.67  ±0.42 

mas/year 

U>y 

=  (AA)tsine  = 

+0.22  +  0.09 

mas/year 

(6) 

(Jjz 

=  -(AA)t  cose  +  (A a)t  = 

-4.57  +  1.29 

mas/year 

where  {A)t  means  the  derivative  of  A  with  respect  to  time.  These  equations 
strongly  suggest  that  the  differences  between  equations  (1)  and  (2)  are  largely 
due  to  the  remaining  rotation  of  the  Hipparcos  reference  frame. 

The  quantities  Aci  and  A 8  are  mainly  determined  from  the  ecliptic  latitudes 
of  the  Moon.  Therefore  they  can  be  affected  by  the  errors  of  Watts’  charts  es¬ 
pecially  in  the  regions  near  the  poles  of  the  Moon.  The  limb  profiles  of  those 
regions  are  now  being  obtained  from  many  lunar  grazing  occultation  observa¬ 
tions  (Soma,  1999).  By  incorporating  such  profile  data,  it  is  expected  that  the 
accuracy  of  the  solution  would  become  much  higher. 


4.  Conclusion 

The  lunar  occultation  analysis  suggests  much  larger  errors  in  the  rotation  of  the 
Hipparcos  reference  frame  than  those  given  by  the  Hipparcos  team. 
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